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Origami and kirigami have emerged as potential tools for the design of mechanical metamaterials 
whose properties such as curvature, Poisson’s ratio, and existence of metastable states can be tuned 
using purely geometric criteria. A major obstacle to exploiting this property is the scarcity of tools 
to identify and program the flexibility of fold patterns. We exploit a recent connection between 
spring networks and quantum topological states to design origami with localized folding motions at 
boundaries and study them both experimentally and theoretically. These folding motions exist due 
to an underlying topological invariant rather than a local imbalance between constraints and degrees 
of freedom. We give a simple example of a quasi-ID folding pattern that realizes such topological 
states. We also demonstrate how to generalize these topological design principles to 2D. A striking 
consequence is that a domain wall between two topologically distinct, mechanically rigid structures 
is deformable even when constraints locally match the degrees of freedom. 


Recent interest in origami mechanisms has been 
purred by advances in fabrication and manufacturing 
[i|-Q, as well as a realization that folded structures can 
form the basis of mechanical metamaterials 0^- The 
ability to identify kinematic mechanisms - allowable fold¬ 
ing motions of a crease pattern - is critical to the use of 
origami to design new deployable structures and mechan¬ 
ical metamaterials. For example, the mechanism in the 
celebrated Miura ori that allows it to furl and unfurl in 
a single motion is also the primary determinant 

of the fold pattern’s negative Poisson ratio UM- Identi¬ 
fying these mechanisms becomes more challenging when 
the number of apparent constraints matches the number 
of degrees of freedom (DOF). 

When there is an exact balance between DOF and 
constraints in a periodic structure, the structure is 
marginally rigid [^, [^. In such a case, new mechanical 
properties such as nonlinear response to small perturba¬ 
tions emerge [iM3- A recent realization is that the 
flexibility of such solids may be influenced by nontrivial 
topology in the phonon band structure [l3:[l8|. Here, we 
show how to extend these topological ideas to origami 
and kirigami. We show that periodically-folded sheets 
may exhibit distinct mechanical “phases” characterized 
by a topological invariant called the topological polariza¬ 
tion, recently introduced by Kane and Lubensky 0 us¬ 
ing a mapping of mechanically marginal structures to 
topological insulators El- The importance of this invari¬ 
ant has er nerg ed in the study of the soft modes of spring 
networks [l8|, and the nonlinear mechanics of linkages 
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[20[ and buckling [2l|. As in these examples, the phases 
in our origami and kirigami structures exhibit localized 
vibrational modes on certain boundaries, and transitions 
from between topological phases are characterized by the 
appearance of bulk modes that cost zero energy. These 
are the hallmarks of topologically protected behavior in 
classical mechanical systems |22l-[^. Topology provides 
a new knob to tune how materials and, as we show here, 
origami and kirigami structures, respond to external per¬ 
turbations. 

We denote by origami, mechanical structures consist¬ 
ing of rigid flat polygonal plates joined by hinges. We will 
first discuss origami with no missing plates or “holes”, 
and then generalize to kirigami, defined to be origami 
where such holes are allowed. We will consider the me¬ 
chanics of origami in the geometrical limit - folds will 
cost zero energy and faces do not stretch or bend. 

To demonstrate the power of our approach, we intro¬ 
duce an example of a ID strip of origami analo gou s to 
the Su-Schrieffer-Heeger polyacetylene model EOT. It 
admits localized modes and stresses determined and pro¬ 
tected by topology, which we realize and characterize in 
experiments. Additionally we show how to generalize 
this to 2D periodic origami sheets, where we have ob¬ 
served a striking property that causes origami without 
holes to have zero topological polarization. We give ex¬ 
amples of hinged structures with holes (kirigami) that 
do admit distinct polarizations and thus can be used as 
building blocks for metasheets with programmable local 
flexibility. 

Quasi-ID origami strip .— We start with a simple 
quasi-ID origami structure. Consider an origami strip of 
zig-zagging rigid quadrilateral plates, depicted in Fig. [TJ 
consisting of a periodically-repeating unit cell of two four¬ 
fold vertices. Each vertex in a cell (labeled W n = 1,2) 


has four creases (Fig. [H^a)), and one DOF |31I that we 
track with the dihedral angles of the bolded crease, /n(i), 
where j indexes the unit cell (Fig. [D^b)). Each adjacent 
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FIG. 1: A quasi-ID origami strip, (a) A unit cell of the fold 
pattern corresponding to the origami mechanism with planar 
angles labeled. Red (blue) creases are mountain (valley) folds, 
respectively, (b) A 3D depiction of a part of the strip with 
folding angles fi(n) labeled, (c) The phase diagram where 
a = TTjS, 13 = I3i =/32 and 7 = 71 = 72 . The colors indicate 
the phase (blue for right-localized, red for left-localized); the 
contours and intensity of color follow the inverse decay length 
1/1 (see legend). Configurations where folds at a vertex be¬ 
come collinear lie on 7 = /3, and the green points along that 
line were constructed in experiment (along this line, ( 7 ,/^) 
and (— 7 , —/3) are related by a rotation in 3D). 


pair of dihedral angles is coupled by the kinematics of the 
intervening vertex. As each vertex contributes a DOF 
and a constraint, this origami structure has marginal 
rigidity. 

We analyze the mechanical response of the origami 
strip by determining its configurations analytically as 
functions of the fold pattern angles 7i and 72 (de¬ 

fined in Fig. [IKa)). We define a generalized displacement 
u{j) = cos/2(j) + 1 . The function u{j) encodes the di¬ 
hedral angle /2 of the right-most fold of unit cell j, and 
satisfies 


u{3 + 1) = K(a,/3i,^2,7i,72)w(i)) (1) 

where 


sm{a — Pi) sin((a — 71) 


sm{a — P2) sin((a — 72) 

_sin(a + j 3 i) sin(a + 71) _ 


_ sin(a + ^^2) sin((a + 72) _ 


( 2 ) 


The derivation is an application of the spherical law of 
cosines and is given in the SI. The fact that u{j) de¬ 


termines u{j + 1) implies that the strip has one degree 
of freedom, globally. Eq. (HD is solved by an exponen¬ 
tial u{j) = 7i(0) exp[j ln(/i:)] with deformation localized 
to one side or the other, following the sign of the inverse 
decay length = In/^. 

The mechanical “phase diagram” in Fig.[T](c) shows the 
values of for patterns with 71 =72 = 7, / 3 i = P2 = P- 
There are two “phases” distinguished by the sign of In 
which is determined here by the sign of 7 + a quantity 
not obviously related to any symmetry-breaking. When 
7 + /3 > (<)0, K < (>) 1 , and by Eq. ([T]), the mechani¬ 
cal response is localized to the left (right) of the origami 
strip. A special role is played by fold patterns with /^ = 1, 
where the decay length diverges and u{j) neither grows 
nor shrinks (denoted by the dashed line). This is pre¬ 
cisely the condition for which a global kinematic mech¬ 
anism exists and the fold pattern is deployable [ 32 |. As 
an example, when 7 = /3 = 0, the strip realizes a row 
of the Miura ori fold pattern, which has a single collapse 
motion. More generically, however, as long as In k never 
changes sign, the deformation in a strip, u{j), is localized 
even if the values of a, vary due to disorder or im¬ 

perfections, i.e. as long as the material remains within the 
same phase. The existence of phases of robust, boundary- 
localized zero-energy deformations separated by critical 
configurations with bulk zero modes suggests that the 
origami strip has topologieally proteeted properties. 

To make the topology explicit, we calculate a topo¬ 
logical invariant of the above phases. Unlike in periodic 
spring networks with marginal rigidity [ 13 , S, Q ^ 
ear analysis is inadequate to capture the topology of the 
origami strip. Coplanar hinges in the flat state are re¬ 
dundant constraints, and this results in extra zero modes 
at linear order which do not extend to higher order. In 
the SI, we derive a rigidity matrix capturing the seeond- 
order deformations of this structure and show that it has 
the same pattern of entries as the Hamiltonian of the Su- 
Schrieffer-Heeger chain of Refs. [ 13,113 Therefore, phases 
of the origami strip are characterized by their topological 
polarization Pt [Slil, defined as a winding number of 
the determinant of the rigidity matrix [ 38 |. Indeed, the 
sign of InK is precisely correlated with the topological 
polarization, and thus the fact that different edges are 
soft or stiff in different phases is a manifestation of the 
bulk-boundary correspondence [HI in this system. While 
topological modes in ID linkages have been found to lead 
to propagating domain walls [HI, IHl , this is not possible 
for our ID strip. In Eq. (j 2 j), k depends only on the fold 
pattern angles not the dihedral angles fj - this 

means that the topological polarization of the unit cell 
cannot change via the zero-energy deformations, which 
would be necessary for propagation. 

To test the consequences of Eq. (HD away from the ideal 
limit, in a structure where faces can bend and hinges can 
twist, a mylar sheet (200 /im thick) is perforated by a 
laser cutter into the desired crease pattern, rendering it 
foldable along lines of perforations. We strengthen the 
facets by sandwiching the mylar sheet between pairs of 
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1 mm thick, plastic plates made of polylactic acid (PLA) 
on a 3 D printer. To mount the plastic plates onto the 
mylar sheet, we use a pushed-in clip design: one facet 
has clips and the corresponding facet has holes. Equiv¬ 
alent holes are cut on the mylar sheet so that the clips 
can be pushed through to meet the holes on the plate on 
the other side of the mylar. An example of the assembled 
origami structure is shown in Fig. Efa). Here, we fixed 
the angle a = 7r/3 and varied 7 = /di = 7 i = /d2 = 72 
to explore the localization of the deformation within one 
phase (with k < 1) [ 39 |. A video camera captured the de¬ 
formation of the strip from above as it was symmetrically 
compressed. The position of each vertex was obtained via 
image analysis, and fit with a 3 D model to reconstruct 
the complete morphology of the origami strip, as shown 
in Fig.[ 2 fb). Finally, the folding angles along the interior 
creases were extracted from the 3 D shape and were used 
to compute the generalized strain u. Fig. [ 2 fc) shows the 
strain as a function of distance along the strip for samples 
with different values of the pattern parameter 7. Observe 
that there is a “soft” edge (cell index 0), where the de¬ 
formation is high, and on the other end a “stiff” edge, 
with low deformation. 

As shown by a semi-log fit (dashed lines in Fig. EJc)), 
the strains decay exponentially at small distances from 
the soft edge. For small 7, the folding angles level off 
to a roughly constant value at larger distances, which 
violates Eq. O- The constant folding angle background 
corresponds to the activation of a mode with uniform de¬ 
formation. This mode is easy to excite as it is the zero 
energy mode at 7 = 0 and thus remains very low en¬ 
ergy for small 7. A deviation from the ideal geometrical 
limit is possible due to the finite flexibility of the facets 
and the finite crease thicknesses. Despite the non-ideality 
of the experimental origami strip, the decay lengths ex¬ 
tracted from the fit are in good agreement with 1/1 = Irk 
(Fig. [ 2 td)), confirming the robustness of our topological 
design principle. 

Two-dimensional origami. — Having established that 
marginally rigid ID periodic origami can exhibit topolog¬ 
ical phases, we now ask whether marginality also leads to 
similar phases in 2D origami. We first characterize the 
class of marginally rigid 2 D periodic origami and show 
that they must have a triangulated crease pattern. To 
avoid trigonometric complexity inherent to a folding an¬ 
gle representation, we model the kinematics of triangu¬ 
lated origami as a central-force spring network with ver¬ 
tices as joints and hinges as springs. Triangles in such a 
network automatically enforce the no-bending constraint 
on the facets. Arbitrary origami can be modeled with 
spring networks, but nontriangular faces require addi¬ 
tional internal springs to remain rigid. 

In this framework, each joint has 3 degrees of freedom 
and each spring adds one constraint, so marginal struc¬ 
tures satisfy E = 3 V where E is the number of bonds 
and V is the number of joints. In a triangulated surface 
without a boundary, each of the E faces is a triangle, 
so 3 E = 2 E. The Euler characteristic y is defined as 



FIG. 2: (a) Localized deformations in an experimental re¬ 
alization of the origami strip (a = 7r/3,7i — ^2 — jdi — 
p 2 = 0.062). (b) 3D reconstruction of the configuration of 
the strip from a flat image (7 = 0.124). (c) Normalized gen¬ 
eralized strain u/uq as a function of distance from the de¬ 
formable boundary (measured in number of unit cells) from 
experiments (shaded curves) with fits to an exponential de¬ 
cay (dashed lines). Each curve shows the average of data 
from 6 — 30 experimental images and has a width equal to 
the standard error. Folding angles /2(1) (related to uq via 
2/0 = 1 + cos/ 2 (l)) at cell 1 varied from 1.06 to 1.45. (d) 
Inverse decay lengths (1//) versus 7 , where data points are 
averages over the fitting coefficients of all images for each 7 
and error bars show 95% confidence bounds. The analytical 
result for = In/^ (Eq.[2|) is plotted as a solid red line. The 
deviations for small 7 arise from a “uniform bending mode” 
(see text). 


X = F — F + F; thus we obtain E = 3 {V — y)- 

Periodic origami structures in 2D have the topology of 
the torus and thus y = 0, which shows that triangula¬ 
tions are marginally rigid. While achieving marginal¬ 
ity in granular packings and glassy networks requires 
some fine-tuning in pressure or coordination, the anal¬ 
ogous origami triangulations arise naturally. Any non¬ 
triangular plate in an origami pattern can be triangulated 
by adding diagonals, and the bending of non-triangular 
plates in real o^ami can be modeled as the addition of 
new creases 0, Q . 

One might now expect a variety of topological phases 
upon changing the angles and lengths of a triangulated 
crease pattern, by analogy with the ID strip. Surpris¬ 
ingly, our calculations indicate otherwise. As discussed 
above, an analysis of the rigidity of flat origami must 
go beyond linear order. To bypass this complication, we 
consider periodic triangulated origami where we break 
the flat-state degeneracy by introducing small vertical 
displacements to the vertices. The linear rigidity and 
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kirigami pattern 



FIG. 3: Topologically protected zero mode (red) in a kirigami 
heterostructure (left green with topological polarization Pt = 
(1,0) / right blue with Pt = (0,0)). Numerically the mode 
depicted has energy nearly indistinguishable from the trans¬ 
lation modes. This is a close-up of a larger periodic 50 x 5 
system, divided into two 25 x 5 domains (two copies in the 
shorter direction are shown). The magnifying glass insets 
show the fine structure of four unit cells of each type, and 
between them is a schematic showing how the quadrilateral 
plates, strips of triangles, and quadrilateral holes are joined 
by hinges. The schematic shows four unit cells, with the lower 
left cell highlighted in green. 


topological properties of such a structure can then be 
expressed in terms of the (Fourier-transformed) rigidity 
matrix R for its associated spring network [l3,[i8|- How¬ 
ever, for all triangulated periodic fold patterns we have 
considered, the function detR(q), a priori a complex¬ 
valued function, is in fact real-valued for all q in the 
Brillouin zone ! Though a proof of this statement for 
all triangulated origami eludes us, extensive numerical 
tests on a large number of distinct fold patterns bear out 
this conjecture. We give details and partial results in the 
SI. 

A consequence of the “reality” property is that the 
winding numbers of Arg detR(q) along any closed 
curves in the Brillouin zone must be zero (when defined), 
and hence the topological polarization Pt must vanish. 
Localized boundary modes for such origami still exist, 
but must be isotropically distributed. Even if the hinges 
in a unit cell break left-right symmetry, the number of 
boundary modes per unit cell on each edge of a finite 
patch is left-right and up-down symmetric. If in fact all 
triangulated periodic origami structures have this prop¬ 
erty, the only way to get an imbalance in the number 
of zero modes at the boundary of origami is by locally 
adding or removing constraints. This behavior contrasts 


with that of the ID strip of origami studied above as well 
as 3D periodic networks and marginal spring networks 
confined to 2D. 

Topological kirigami .— Thus the question remains: do 
there even exist 2D periodic hinged structures with a 
nonzero topological polarization? The answer is yes, but 
we must go beyond origami to kirigami^ folded structures 
with holes. There is a simple way to generate marginal 
kirigami from triangulated origami. Cutting out an adja¬ 
cent pair of triangles removes one bond from the associ¬ 
ated spring network, eliminating a constraint. Likewise, 
merging two triangles into a rigid quadrilateral plate adds 
a constraint. We therefore modify a triangular lattice by 
cutting and merging twice, resulting in a structure with 
two quadrilateral plates and two quadrilateral holes per 
unit cell (top center of Fig. [3]). Now det R(q) is complex¬ 
valued, and by randomly perturbing a flat realization, 
we find the “green” (left) and “blue” (right) structures 
depicted in Fig. O which have Pt = (1,0) and (0,0), re¬ 
spectively (see SI for more details). With free boundary 
conditions, the boundary soft modes in the green kirigami 
are polarized to the Px edge (analogous to the ID strip 
and in contrast to the blue kirigami and all triangulated 
origami structures we tested). 

Finding the green kirigami answers the question above 
positively, and we leave a determination of the possible 
phases that can occur in the modified triangular lattice to 
future work. A full characterization will likely be difficult 
due to the high dimensionality of the realization space 
(c.f. Ref. [ 3 ^ which shows the complexity of the phase di¬ 
agram in a simpler mechanical system). We thus switch 
gears and present an example of localized modes designed 
into a kirigami “heterostructure” to illustrate the power 
of our techniques. In Fig. [3l we show a kirigami struc¬ 
ture that exhibits zero modes localized at a domain wall 
(one per unit cell) between the two kirigami structures 
described above. The zero modes render the heterostruc¬ 
ture flexible in the vicinity of the domain wall (the mode 
depicted leads to out-of-plane wrinkling), while keeping 
it rigid further away. By contrast, a domain wall be¬ 
tween distinct patterns with equal polarization has no 
such localized modes (see SI). In general, domain walls 
between structures with different polarizations create ei¬ 
ther “soft” lines along which the system easily deforms, or 
“stressed” lines which first bear the loads under applied 
strains [2^. Similar effects may arise at point defects in 
otherwise uniform polarized structures . 

Outlook .— We have demonstrated that origami and 
kirigami structures are characterized by a topological po¬ 
larization that classifies the ways that a marginaiiy rigid 
foid pattern can be floppy ciose to its boundaries. Our 
resuits give strong constraints on the types of boundary 
modes that can be created in origami and wiii guide the 
design of foid patterns that achieve a targeted mechan- 
icai response. In the design space of geometric reaiiza- 
tions, two structures with different poiarizations must 
be separated by giobaiiy flexibie, i.e. depioyabie, reaiiza- 
tions. Thus not oniy can structures with distinct phases 
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be combined in real space to form domain walls with 
useful functionality, but also they can be used to find 
deployable patterns in design space. These realization 
spaces are high dimensional in general, so the problem 
of determining simple rules to create a given polarization 
remains open. 
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FIG. SI-1: Spherical polygons for a four-fold vertex. Angle 
variables are different from those used in the main text. 


Supplementary Information 

Appendix A: Vertex kinematics and spherical 
polygons 

The intersection between an origami vertex and a 
sphere traces out a spherical polygon on the sphere sur¬ 
face [sil (Fig. ISI-ip . For a four-fold vertex, this poly¬ 
gon is a quadrilateral whose side lengths are inherited 
from the angles between adjacent folds in the fold pat¬ 
tern (the unlabeled fourth angle in Fig. ISI-ll is 7 “ = 
27V — — P~ — 7 ^)- The relationship between dihedral 

angles is 

coscos 7 ^ + sinsin 7 ^ cos/^ = (SI) 
cos P~ cos 7 “ + sin P~ sin 7 “ cos f ~. 

Since /3+ + 7 + = 27r — — 7 “, Eq. ([ST]) becomes 

sin sin 7 ^ (cos /^ + l) = sin P~ sin 7 “ (cos f~ + l)(S 2 ) 

A similar relationship between 0^ and 0~ can be de¬ 
rived as 


sin 7 + sin 7 (cos -|- 1 ) = sin P^ sin P (cos 0 + IX-SS) 


The relationship between /+ and however, is more 
complex fFig. ISI- 2 ]) . Define e to be the angle between 0^ 
and 0~ given by 

cos e = cos P^ cos 0^ + sin p^ sin 0^ cos /+. (S4) 


and such that e = P^ + 7 ^ when the vertex is flat. Then 


(9+ = 


cos 


± cos 


f cos P^ — cos e cos 7 + 


sm e sm 7 ^ 

^ f co^P~ — cose cos 7 “ 
sin e sin 7 “ 


(S5) 



FIG. SI- 2 : The angle 0^ versus for two cases: (solid) = 
P~ = 7 r/ 3 , 7 + = 7 “ = 27r/3 and (dashed) P~^ = 7 r /3 — tt/IO, 
P~ = 7 r /3 + tt/IO, 7 ^ = 7 ~ = 27r/3. In each case, there are 
two distinct branches joined only at the flat state. Note that 
one of the branches for the solid curves lies along /+ = 0 . 


Either sign in Eq. ([S5]) satisfies the geometric con¬ 
straint equations; this ambiguity arises from the differ¬ 
ent branches of cos“^. Yet, the branches are physical, 
and correspond to different ways of choosing which folds 
are “mountain” and which folds are “valley” folds. In 
particular, there is a branch in which /+ and f~ are 
both either mountain or valW folds and one in which 
they have the opposite sense Q. In either case, however, 
cos remains the same and the exponential solution for 
u = 1 cos / in Eqs. 1 and 2 of the main text is valid. 

One special case bears additional analysis. When the 
vertex has 7 ^+/^^ = tt, the folds 0^ and 0~ are collinear. 
In that case, one branch of the vertex has both /+ and 
f~ as mountain (valley) folds. Eor the other branch, 
however, /+ = /“= tt is strictly satisfied and = 
0~. That is simply the case that a vertex is folded along 
the line shared by and 0~ ^ and this prevents the two 
adjacent folds from folding at all. This motivates the 
choice of variables discussed in the main text, for which 
P^ = P and 7 ^ = 27r — a ^ 7 . Then p = ^ precisely 
when 0^ are collinear. 

Appendix B: Second-order rigidity matrix for ID 
strip 

We label each unit cell with an index j and the two ver¬ 
tices within a unit cell with n = 1, 2. Eor convenience, 
we also express the dihedral angles in terms of their dif¬ 
ference from TT, Sf^{j) = fnU) ~ tt for the folds of 
vertex n within unit cell j. Each vertex has one DOE, 
which we will parameterize with S0^{j). Since there is 
a two vertex unit cell, we combine these parameters into 
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two vectors, 


and 


f(i) = 


/r(i) \ 

fiU) 

/2 U ) 

ftU)) J 



(S6) 


(S7) 


Using spherical trigonometry, we compute = Js(jf), 
where 


f Ai ^ (JiBi 0 ^ 

0 

0 A 2 + <^2-^2 ^ 

0 A+ + £72^2+ / 


(S8) 


(Jn = ±1 is an arbitrary choice of branch in configuration 
space the vertex should be folded. 


± ^ sin(Q =F 7n) 
" sin(2a) 


(S9) 


and 


_ 1 / sin(Q: - 7n) sin(Q: + 7^) sin(Q: =F ^n) 

" sin(2Q;) y sin(Q; ± 

(SIO) 

For this paper, we have always assumed that (7^ = 1. 

To complete our description, we must enforce the lin¬ 
ear constraints Sfi{j + 1) — = 0 and Sf^ij) — 

= 0. In Fourier space, these constraints lead to 
an equation of the form 


R(g)s(^) = 0. (Sll) 


for the infinitesimal deformations of the origami strip, 
where 


R(^) = 


{A-^B-)e^^ -{A^^B})\ 

^2 + ^2 ^ J 


(S12) 


Appendix C: Real-valuedness of the determinant: 
a hidden symmetry 

We will restate the observed “reality” property here 
more formally. For all periodic triangulations we have 
checked, there exists a choice of unit cell so that the de¬ 
terminant of the rigidity matrix det R(^ is real valued for 
all ^in the Brillouin zone. The consequence of this is that 
with such a choice of unit cell, the phase of det R(g) can¬ 
not wind and so the topological polarization must vanish. 
There are also obvious consequences for the set of bulk 
zero modes which we will not discuss here. In Fig. ISI-31 
we display a sample calculation and visualization of this 
propety for a small periodic triangulated origami, repre¬ 
sented by a spring network. 

This real-valuedness can be expressed in the following 
equivalent way. First, define and ^2 = ^ 



(a) (b) 


-^ (-17085663 + 33649854^1 - 16524255^? + 1545166^2 

37647680002;i2;2 

- 1431150^122 - 88999821^22 + 6162382^22 + 108292 | 

+ 29283730212| - 583495022?2^ + 292837302?2| 

+ 108292?2| + 616238212| - 8899982?2^ - 14311502?2| 
+15451662?2f - 165242552?2| + 336498542?2^ - 170856632?2|) 


(c) 




Zi 

4 

Zi 

4 

^2 

0 

0 

-16524255 

33649854 

-17085663 


0 

616238 

-889998 

-1431150 

1545166 


10829 

29283730 

-58349502 

29283730 

10829 

^2 

1545166 

-1431150 

-889998 

616238 

0 


-17085663 

33649854 

-16524255 

0 

0 


(d) 


• • • 

1 .... 

0 • • • • • 

-1- • 

-2 • 0 • 

-2-1012 

(e) 


FIG. SI-3: Illustration of the reality property for triangulated 
origami, (a) A periodic triangulation based on the triangular 
lattice with four vertices, 12 edges and 8 triangles. The vertex 
coordinates are (0,0,0), (0,1,1/5), (1,1/7, 0), (3/4,1,—1/7) 
and the lattice primitive vectors are (2,0,0) and (0,2,0). 
Above is a side view and below is a top view. Four unit cells 
are shown, with one highlighted in green, (b) The Fourier 
transformed rigidity matrix Il(zi,Z 2 ), where zi = and 

Z 2 = Every row corresponds to an edge in the trian¬ 

gulation. (c) The Laurent polynomial detR(zi,Z2). (d) The 
table of coefficients of det R(zi, ^2), organized by powers of 
and Z 2 . The number 3/3764768000 has been factored out of 
all entries, (e) A graphical representation of the coefficients 
of det R(zi, 2:2). Red dots are positive coefficients, blue dots 
are negative, and the area of the disks is propoertional to the 
absolute value of the coefficients (with a minimum disk size 
imposed for visibility). 


and observe that det R(2:i, 2 : 2 ) is a Laurent polynomial in 
those variables. Next, note that det R(^ is real if and 
only if: 

det R(^ = det R(^ 
det R(2 :i , 2 : 2 ) = det R( 2 :i, 2 : 2 ) 

= det R(zi, ^ 2 ) 

= det ^ 2 ^)' 

In the third line, we use the fact that the coefficients of 
this polynomial are real, and in the final line, we use the 
fact that all q in the Brillouin zone, zi and 2:2 lie on the 
unit circle and so Zi = z^^. This shows that for all 
the coefficients of the zfz^ and z^^z^^ terms must be 
equal in order for this polynomial to be real-valued on 
the unit circle. 
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A change of unit cell merely changes det R by a factor 
Zi^Z 2 where m and n are integers related to the change 
in the number of bonds crossing unit cell boundaries [13 • 
If one imagines the coefficients of detR as a table (Fig. 
ISI-Sr d)). then such a factor does not change the values 
in the table, but merely shifts the positions of the entries 
in the table by a vector (m,n). Therefore, the “reality 
property” is equivalent to the polynomial det R(2;i, 2 ^ 2 ) 
being “2D palindromic”, i.e. 

det R(zf Z 2 ^) = Z 1 Z 2 ^ 2 "^)’ (S13) 

for some integers k^l. Equivalently, the table of co¬ 
efficients of det R(2 ;i, 2 ^ 2 ) is symmetric under inversion 
through some point {kl2Jl2). 

Examples and tests performed: We generated small trian¬ 
gulations of a torus via flipping edges from the triangular 
lattice and in all cases checked found that the the polyno¬ 
mial det R(2 :i, 2 ^ 2 ) was a palindrome. To ensure that loss 
of precision from floating-point operations were not an 
issue, we did these computations in rational arithmetic 
in the Mathematica computer algebra system. 

Specifically, we began with periodic triangulation 
“seeds” displayed in Fig. ISI-41 and performed random 
edge flips (transforming two adjacent triangles into two 
different triangles on the same four vertices) repeatedly 
onto randomly chosen edges. This procedure is illus¬ 
trated in Fig. ISI-4r ab 

For every such crease pattern generated, we used ver¬ 
tex coordinates made by randomly perturbing a fiat em¬ 
bedding of the seed triangulation. These were rounded 
to the nearest rational number with denominator lower 
than 10. The Laurent polynomial det R(2;i, 2 ^ 2 ) was then 
calculated exactly and the coefficients of different terms 
were compared to check whether the polynomial was in¬ 
deed palindromic as in Eq. (jS13p . In some cases, a flip 
would cause this polynomial to vanish (due e.g. to a de¬ 
generate edge connecting a vertex to itself). Afterwards, 
it is quite likely that further flips would produce trian¬ 
gulations with zero determinant. While the zero polyno¬ 
mial is of course real-valued, if this occurred, in order to 
get more interesting tests, we rejected this triangulation 
and flipped another edge instead. Regardless, for all tri¬ 
angulated origami checked we found that Eq. (1^ was 
always satisfied. In Table [I] we show a summary of the 
tests performed. 

Partial result: We can show that subdividing a trian¬ 
gle in triangulated origami preserves the property that 
det R(^ remains real. To see this, note that such a sub¬ 
division creates one new vertex and three edges. With 
an appropriate choice of unit cell we can ensure that this 
vertex is not on the boundary and thus does not pick up 
factors of or 

The new vertex introduces three columns into R(^ 
whose only nonzero entries can be in the rows corre¬ 
sponding to the three new edges. By reordering the ver¬ 
tices and columns, one finds that the determinant takes 
a block upper-triangular form with a 3 x 3 block on the 
diagonal from the new additions (which we call Rn(^) 



6x4 seed 


(b) 


/ 



/ 


1 


7 

/ 


/ 



7 

7 

7 

/ 

/ 




7 


7 

/ 

z 


7 




7 

y 



7 


7 

7 

7 



/ 



7 

7 


7 

7 


one edge flip 

(a) 




/ 




7 

7 


7 


^7 

/ 



7 

7 






7/ 




/ 

7 







7/ 



/ 

/ 







7 


y 


/ 







7 


7^ 

y 

/ 

/ 


7 




7 


7 

7/ 

y 

/ 


/ 


7 


7 


7, 

7 

7/ 

^7 

/ 

/ 

7 

7 

7 

/ 



7 



y 

/ 

/ 

7 

7 

7 



7 



7. 

y 

/ 

/ 

7 

7 

7 


7 




7. 

y 


/ 

7 

7 

7 

7 




7 

7/ 

y 


/ 

7 

/ 

7 




7 

7 

7/ 



6x6 seec 


(C) 



(d) 


FIG. SI-4: Seed triangulations and edge flips. In each im¬ 
age, four unit cells are shown, with one highlighted in green. 
The triangulations are only shown with planar embeddings for 
clarity, the actual tested triangulations were deformed into 3D 
with random perturbations, (a) left: The 4 by 4 seed trian¬ 
gulation with one highlighted edge, center: the triangulation 
resulting from the edge flip of the red edge, right: a triangula¬ 
tion resulting from 29 random edge flips on the seed, (b) The 
6 by 4 seed triangulation, (c) The 6 by 6 seed triangulation, 
(d) The 8 by 8 seed triangulation. 


seed |F| \E\ \T\ (runs, flips) 

4 by 4 16 48 32 (100, 1000) 

6 by 4 24 72 48 (50, 500) 

6 by 6 36 108 72 (20, 200) 

8 by 8 64 192 128 (10, 100) 

TABLE I: Summary of tests performed of the “reality” prop¬ 
erty on different families of periodic triangulated origami. 
Each line of the table corresponds to a different seed (de¬ 
picted in Fig. ISI-4D . and displays the number of vertices |F|, 
edges \E\ and triangles |T| in the seed and shows the number 
of runs started as well as the number of random flip opera¬ 
tions performed. Note that a flip preserves |F|, \E\ and |T|, 
and that the computations require computing determinants 
of symbolic \E\ x 3|F| matrices. In all runs, every single tri¬ 
angulation generated had the reality property, as described in 
the text. 


as well as Ro(^ as another block on the diagonal. By 
a well-known property of the determinant, this implies 
that detR(^ = detRn(^ detRo(^. 

If we were able to prove an analogous result for edge 
flips, then we would have shown that the palindromic 
property holds for all triangulations. 

Appendix D: On the design of topological 
kirigami 

We found topologically polarized kirigami via the fol¬ 
lowing steps. We first constructed the pattern of plates 
and hinges as follows. As in Appendix C, we began with 
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kirigami I 


kirigami II 





Pt = (0,0) 


FIG. SI-5: Schematics of some structures created while de¬ 
signing the kirigami. In each image, four unit cells are shown, 
with one highlighted in green. These images only show the 
pattern of hinges and plates, and do not represent 3D em¬ 
beddings. (a) A triangulated origami based on the triangular 
lattice. Two families of creases are highlighted in red and 
blue which are moved to create the next structure, (b) The 
checkerboard lattice of tetrahedra. This structure consists of 
tetrahedra joined by universal joints at their vertices, (c) The 
final kirigami structure. Created by placing rows and columns 
of additional triangles between tetrahedra, and then replacing 
the tetrahedra with flat quadrilateral plates. 


a periodic spring network with the combinatorics of the 
standard triangular lattice (drawn in Fig. ISI-bf a) as a 
square lattice with all NE diagonals added). This pattern 
models a triangulated origami structure, which seems to 
have a real determinant. To try to break this symmetry 
we attempted to add holes in the unit cell. 

To add holes without changing the marginal rigid¬ 
ity property, we removed one NE diagonal bond from 
a square A and placed it in an adjacent (already occu¬ 
pied) square 5 as a NW diagonal bond. This yields a 
hole in A and a new tetrahedron in B. If one does this 
for all squares in a checkerboard pattern, one ends up 
with the so-called checkerboard lattice of tetrahedra (Eig. 
ISI-5lf b)). This is the simplest 2D periodic structure in 3D 
which admits geometric realizations with distinct topo¬ 
logical polarizations. However, the tetrahedra are joined 
at vertices, not along hinges, so it is neither origami nor 
kirigami. 

To remedy this, we may add rows and columns of trian¬ 
gulated strips between the holes and the tetrahedra. Now 
we do have a structure made of solid bodies joined along 


FIG. SI-6: Three kirigami unit cells with different topologi¬ 
cal polarizations. A top view and side view of the 3D unit 
cells are shown above exploded versions that have been cut 
into strips and then flattened into the plane. The dashed lines 
indicate how the strips are to be joined to create the 3D struc¬ 
ture shown above and the faces labeled by letters indicate the 
periodicity of the structure. 


hinges, but the bodies are tetrahedra and not yet flat 
plates. Flat plates can be modeled by a spring network 
with pyramids. By adding an additional vertex above 
the plate and attaching it to the four vertices of the we 
create a minimally rigid body which is attached to the 
structure by four coplanar hinges. The final pattern of 
hinges and plates is shown in Fig. ISI-bf c) (pyramids not 
shown for clarity). 

Given this design, we still have to find a suitable re¬ 
alization of this structure as a 3D object. We did this 
by beginning with a fiat embedding of a unit cell and 
perturbing the vertices randomly, subject to the condi¬ 
tion that the quadrilteral plates remain planar. These 
randomly generated unit cells were then checked to de¬ 
termine if they had a topological polarization. Fig. ISI-61 
shows the unit cells of three kirigami structures with this 
crease pattern but with differing topological polariza¬ 
tions: kirigami I has polaarization (1,0), and kirigami 
II and III have polarization (0,0). Their 3D structures 
are also displayed in an “exploded” view, where they have 
been cut apart into flattened strips which must be bent 
in spance and then glued together to form the kirigami 
structures. Unit cells I and II were used to construct the 
domain wall in Fig. 3 in the main text. 

In Fig. I SI- 71 we show kirigami structures consisting of 
pairs of unit cells from Fig. ISI-61 The domain walls do 
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not contain any extra or missing bonds, so there are 
no differences in the local counting of degrees of free¬ 
dom. These kirigami structures are rectangular sheets 
with periodic boundary conditions in both directions, so 
one might expect the lowest energy vibrational modes to 
be bulk bending modes. This is indeed the case when 
the two unit cells used have the same polarization vec¬ 
tor (Fig. ISI-7r a)). However, when a structure is made 


of unit cells with two distinct polarizations, one of the 
domain walls carries localized very low-energy modes 
(Fig. ISI-7r b) and Fig. 3 in the main text). The eigen¬ 
value corresponding to the energy of the topological mode 
is uj‘^ = 1.42 X 10“^^ (comparable to the eigenvalues cor¬ 
responding to the translational zero modes), whereas the 
eigenvalue corresponding to the bending mode is much 
larger, = 1.82 x 10“^. 
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FIG. SI-7: Composite kirigami systems constructed from unit cells I, II and III from Fig. lSI-Gl and their lowest energy eigenmodes. 
Both structures depicted consist of 50 unit cells (along the x-direction) by 5 unit cells (along the ^-direction), which are divided 
into two domains each of size 25 by 5. Two copies in the shorter direction are shown in the upper figures, and the displacements 
in the mode are overlaid as red lines. The structures are periodic in both directions, so there are domain walls at x = 0 = 200 
and X = 100. The domain wall at x = 100 is marked with dashed lines. The lower plots show the x-, y- and ^-components (in 
blue, green, and red, respectively) of the displacement vectors of all vertices as a function of the x-positions. (a) Composite 
structure of kirigami II and III, both of which have polarization (0,0). The lowest eigenmode has eigenvalue = 1.82 x 10“® 
and resembles a bulk bending mode, (b) Composite structure of kirigami I (polarization (1,0)) and II (polarization (0,0)). The 
lowest eigenmode has eigenvalue = 1.42 x 10“^® and is highly localized near the domain wall at x — 100. 













